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Abstract 

Predictions for electron induced proton knockout from the pi/2 and P3/2 shells 
in ^^O are presented using various approximations for the relativistic nucle- 
onic current. Results for the differential cross section, transverse-longitudinal 
response (Rtl) and left-right asymmetry Atl are compared at \Q'^\ = 0.8 
(GeV/c)^ corresponding to TJNAF experiment 89-003. We show that there 
are important dynamical and kinematical relativistic effects which can be 
tested by experiment. 
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Under quasielastic conditions the (e, e'p) reaction can be treated with confidence [|T| via 
the impulse aproximation in which the detected proton absorbs the entire momentum {q) 
and energy {uo) lost by the electron. Detecting the outgoing proton in coincidence with 
the scattered electron allows one to determine the missing energy {Em) and momentum 
{pm) in the reaction and thus to provide detailed information on the energies, momentum 
distributions and spectroscopic factors of bound nucleons [|l],H . Until recently low- -Em. data 
were concentrated at pm < 300 MeV/c. Now higher p^-regions are being probed at small 
Em under quasielastic conditions , yielding new information in a regime where two-body 
currents (meson exchange currents and A-isobar contributions) can be safely neglected . 

Most theoretical work on (e, e'p) has been carried out on the basis of nonrelativistic 
approximations to the nucleon current. Specifically, the standard distorted wave impulse 
approximation (DWIA) [|I| has been extensively employed to analyze (e, e'p) data, using 
nonrelativistic current operators with bound and scattered proton wave functions deduced 
from phenomenological nonrelativistic potentials. However, DWIA analyses have met two 
major difficulties: a) The spectroscopic factors extracted from low-p^ data are too small 
compared with theoretical predictions For instance, the extracted occupations of 3si/2 
and 2(^5/2 orbits in ^osp]^ ^^^^ ^ g while theories on short-range correlations [0 predict 
at most a 30% reduction of mean-field occupations for levels just below the Fermi level, b) 
DWIA calculations compatible with the low-p^ data predict much smaller cross sections 
at high-pm than those experimentally observed 0. Although short-range correlations are 
expected to increase the high-momentum components, their effect is negligible ||^ at the 
small missing energies of these high-p^ data, and effects of long-range correlations have 
been invoked 0. 

In recent years the relativistic mean-field approximation has been sucessfuUy used for 
the analyses of both low-p^ @HI1I ^^"^ high-p^ data. In the relativistic distorted-wave 
impulse approximation (RDWIA), the nucleon current 



(1) 

is calculated with relativistic ipB and ipF wave functions for initial bound and final outgoing 
nucleons, respectively. Jp^ is the relativistic nucleon current operator of ccl or cc2 forms 
as in |jl5|. The bound state wave function is a four-spinor with well-defined parity and 
angular momentum quantum numbers and is obtained by solving the Dirac equation with 
scalar-vector (S-V) potentials determined through a Hartree procedure from a relativistic 
Lagrangian with scalar and vector meson terms |]TB|. The wave function for the outgoing 
proton is a solution of the Dirac equation containing S-V global optical potentials ||T^ for 
a nucleon scattered with asymptotic momentum pp. In contrast to DWIA, in RDWIA the 
analyses of individual nuclear shells were done with no fitting parameters other than the 
spectroscopic factors |llOHl3|. The RDWIA spectroscopic factors obtained are larger than 
the DWIA ones |Tl|-p!^ and are valid both for low- and high-p^ data [|1^ 
mentioned 3si/2 and 2^3/2 shells in ^°^Pb, values of — 0.7 have been obtained 
together with reasonable agreement at high-p, 
We have recently studied 



for the above- 
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relativistic treatment of the nucleon current 



T8| , p!9| the effect on the individual response functions of the 
Focussed on the relativistic plane wave im- 



pulse approximation (RPWIA) we showed |jT8| that the TL response is very sensitive to the 
negative-energy components of the relativistic bound nucleon wave function. We have also 
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shown that, for the j = I ±1/2 spin-orbit partners of a given shell, this sensitivity is 
much larger for the j = / — 1/2 than for the j = 1 + 1/2 case. Strong sensitivity of the TL 
response to relativistic corrections was earlier found for d{e, e'p) [ pO[| . 

A certain degree of controversy surrounds the TL response measured in exclusive 
quasielastic electron scattering from the least bound protons in several nuclei (^^C, ^^O, 
^°^Pb): in some cases ||2l|] large deviations from standard DWIA calculations appear, while 
in others the data are close to the calculations. New data on the Rtl response for proton 
knockout from the lpi/2 and lp3/2 orbits of ^^O are expected soon from an experiment car- 
ried out at Jefferson Laboratory (TJNAF) [g^f at \Q'^\ = 0.8 (GeV/c)^ and perpendicular 
kinematics. The purpose of this work is to show for this case that there are important 
kinematical and dynamical relativistic effects. 



For the bound state wave functions we use the parameters of the set NLSH in p3 



We have also obtained results with the older HS set p6ip4| , as well as with the newest 
NL3 one and found very similar results. For the scattered proton wave function, we 
use the energy-dependent, A-independent, potentials derived by Clark et al. [O for ^^O. 



Again, we checked that using other relativistic optical potentials does not alter our results 
and conclusions. The kinematical conditions follow closely those of the TJNAF proposal 
2^: beam energy 2445 MeV, |g| ~ 1 GeV/c, uj ~ 439 MeV, corresponding to quasielastic 



kinematics (a; ~ \Q'^\/2/M), and proton kinetic energy of about 427 MeV; the angle Op of 
the outgoing proton with respect to q is varied and thus the missing momentum varies as 

Pm = W + - 2|Pf||'?1 COs6'i7'. 

Figure 1 shows the differential cross section, Rtl response and TL asymmetry {Atl) 
for (top panels) and P3/2 (bottom panels). Rtl and Atl are obtained from the cross 
sections measured at (pp = 0° and (pp = 180° with the other variables {UjQ"^, Em,Pm) held 
constant, where (pp is the azimuthal angle of the scattered proton (we follow the same 
convention as in p2|). Relativistic calculations using the ccl and cc2 current operators are 
shown by solid and dotted lines, respectively. The Coulomb gauge has been used throughout 
and we have checked that, as in our previous work (see in particular fl^), the Landau 
and Coulomb gauges produce similar results. On the contrary, the Weyl gauge produces 
important deviations and tends to give unrealistic enhancements of the longitudinal current 
ITSl . We use a spectroscopic factor Sa = 0.7 as obtained for ^osp]^ which in RDWIA also 
matches the low-p^ data P,|26[] for the shells discussed here; the spectroscopic factor simply 
scales down the curves for differential cross sections and Rtl while leaving Atl unchanged. 

We can divide the differences between this fully relativistic approach and the standard 
nonrelativistic one into two categories: i) Effects due to the fully relativistic current opera- 
tor, i.e., 4x4 matrix structure of the 4- vector current operator, compared to the 2x2 matrix 
structure of the nonrelativistic current operator which is usually expanded in p/M. We call 
these effects kinematical as they are in principle independent of the dynamics introduced 
by the nuclear interaction, ii) Effects due to the differences between relativistic and non- 
relativistic nucleon wave functions, which of course depend not only on the 4-spinor versus 
2-spinor structure, but also on the potentials used in the respective Dirac and Schrodinger 
equations. Thus we call these effects dynamical. 

Further insight into these differences can be gained recalling the steps taken in nonrel- 
ativistic approaches. First, the one-body current operator is expanded in a basis of free- 
nucleon plane-waves, which amounts to a truncation of the nucleon propagator that ignores 
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negative-energy solutions of the free Dirac equation. This is common to PWIA and DWIA. 
Second, to take into account (spin- dependent) final-state interactions or distortion effects, 
in DWIA the current operator is transformed into a 2 x 2 matrix (usually involving p/M 
expansions) to calculate the nucleon current as the matrix element between the nonrela- 
tivistic bound {xb) and scattered (xf) wave functions. This second step is not needed in 
PWIA, which can then better incorporate the relativistic kinematics, but misses important 
absorption effects. A "relativized" form of the 2x2 current operator was proposed in P7|J2^ 
to optimize the p/M expansion. 

One may then identify two types of relativistic dynamical effects: 

ii — a) Dynamical effects coming from the difference between the upper components of tpp 
{ips) and the solutions xf (xb) of the Schrodinger equation. Assuming equivalent central 
and spin-orbit potentials, this difference stems from the well-known Darwin term. Provided 
that the relativistic dynamics are known, one can deduce the Darwin term and construct 
an equivalent bi-spinor wave function x fo include its effect in the nonrelativistic nucleon 
current, thus removing this source of difference between relativistic and nonrelativistic results 



. This is done for instance in ref . p9| . The influence of this term on (e, e'p) observables 



has been demonstrated in several works 1 11, 12,30, 31 1. It appears to be the main dynamical 



relativistic effect in the cross section in the low-p^ region |TT] , pr2[ , and is important for the 
correct determination of the spectroscopic factor from low-p^ data. Its omission reduces the 
spectroscopic factor by 15-20%. It is included in all calculations presented here. 

ii — b) The other dynamical effect is due to the negative-energy components of the 
relativistic ipB, "ipF wave functions. Starting from Schrodinger-like solutions x may at 
best construct properly normalized four-spinors of the form 

to calculate the relativistic nucleon current. However, this spinor lacks the dynamical en- 
hancement of the lower component of the Dirac solution due to the relativistic S-V potentials. 
This dynamical enhancement is contained in the negative-energy components of the rela- 



tivistic ipB {"^f) solutions and influences {e,e'p) observables in the high-p^ regions |T2|,|T8 
A discussion of how this effect may be incorporated in nonrelativistic formulations based 
on 2 X 2 matrix, p/M expansions of current operators modified to include the effects of the 
relativistic S-V potentials can be found in p2 |. 



The role of the negative-energy components can be seen in Fig. 1 comparing the solid 
with the short-dashed lines. The dashed lines show the results obtained with the ccl current 
operator when the negative-energy components are projected out, i.e., the nucleon current 
is calculated as 

J'Ni++)i^, Q) = j dM^\p + q}>{u;, g)#(p-), (3) 
where ipj^^ {'^Pf'^) is the positive-energy component of tpB (iPf), i-e., 

^^b\p} = ^^MMp), A(+)(p1 = (4) 

with = {E,p) and E = y/p^ + M'^ (similarly for "ipj^^)- The difference between the 
solid and short dashed lines is due to the dynamical enhancement of the lower components 
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which is contained in the current of eq. (||), but not in eq. (^). This effect is more 
visible than that introduced by the theoretical uncertainty due to the choice of ccl (solid 
line) or cc2 (dotted line) operators. It is important to realize that the positive-energy 
projectors inserted in eq. (||) depend on the integration variable p. One may attempt to 
neglect this dependence by using projection operators corresponding to asymptotic values 
of the momenta, i.e. projectors acting on ipp and ipB respectively, with Pp = {Ef.Pf)-, 
Pp — the asymptotic four-momentum of the outgoing and bound nucleon respectively, 



with = {LU,q) and u = Ep — y {pf — o)"^ + M"^- We refer to this approach as asymptotic 
projection (Jas)- The results corresponding to this approximation are shown by long dashed 
lines in Fig. 1. They are obtained with the ccl operator and are very similar for cc2. 

The differential cross sections for \pm\ < 300 MeV/c are similar, but show a substantial 
dependence on the negative-energy components for \pm\ > 300 MeV/c for either the pi/2 or 
P3/2 shells. Note also that the cross sections obtained with positive-energy projected wave 
functions are more symmetrical around Pm = than the RDWIA results. Therefore, the 
effect of removing the negative-energy components shows up more in Rpi and App (see 
middle and right-hand panels of Fig. 1). Particularly interesting is the oscillatory structure 
of the fully relativistic result for Atl- This characteristic is preserved by the positive-energy 
projection method of eq. (|), but not by the method of asymptotic projection. Note that the 
dependence on the dynamical enhancement of the lower components is stronger for the pi/2 
Rtl response than for ^3/2, a feature that was first seen in RPWIA [0 and that persists 
in RDWIA. On the other hand, the asymptotic projection severely modifies Atl for both 
orbitals. We notice that the App calculated with Jas approach are very similar to the ones 
obtained in |2^, as it is similar to the EMA (noSV) in said reference. At low momentum this 
approach lies close to the fully relativistic one and to the Jproy ones, but beyond Pm — 200 
MeV/c it gives completely different results. The oscillating trend of the App calculated in 
RDWIA is confirmed by the preliminary data [^] and agrees qualitatively with previous 
calculations by Van Orden p2| . 

Other relativistic effects can be seen in Fig. 2, where we compare RDWIA results on 
Atl (left panels) and Rtl (right panels) to nonrelativistic approaches at various levels. To 
minimize the differences we have used the cc2 current operator and nonrelativistic scattered 
wave functions obtained from Dirac-equivalent Schrodinger equations [|I^]. This ensures 
that the nonrelativistic wave functions correspond to the upper components of the relativistic 
ones, containing in particular the Darwin term. For the nonrelativistic bound wave functions, 
we used the ones in . In p7| , p8| new approximations to the on-shell relativistic one-body 



current operator were developed to take better account of relativistic kinematic effects in 
nonrelativistic calculations. In particular, the charge density contains a spin-orbit correction 



that affects Rtl 1131 • Fig. 2 we show by long dashed lines the results obtained with the 
"relativized current" and by dotted lines the results obtained with that relativized current 
when the spin-orbit correction term to the charge- density is neglected. One can see that 
the spin-orbit correction has a very large effect on Rtl and Atl- Its omission causes large 
deviations from the relativized current results. Using the DWEEPY |^3[ code we have 



obtained similar results to the dotted lines in Fig. 2. 

The short dashed lines in Fig. 2 are results obtained with the relativistic current and 4- 
spinors constructed as in eq. @j from the nonrelativistic bound and scattered wave functions. 
In this way, the relativistic kinematics are fully taken into account as well as the dynamical 
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effect in the upper component, namely the effect of the Darwin term; only the dynamical 
enhancement of the lower components is missed. This is why these results (short-dashed 
lines) for A^l and Rtl are much closer to the fully relativistic results shown by the solid 
lines. We see also that a large oscillation of Atl can be recovered in the nonrelativistic 
approach. The much smaller oscillation of Atl is a distinctive feature of the Jas results, 
and it seems to be ruled out by the experiment. We also note that, while the effect of the 
dynamical enhancement of the lower components is larger in pi/2 than in shells, the 
effects of relativistic kinematics are of the same order in both shells. 

In conclusion, we have identified two types of relativistic effects on Rtl and Atl- One 
is of kinematical origin, and has a large contribution from the spin-orbit correction to the 
charge density, and other is of dynamical origin. The latter is mainly due to the enhancement 
of the lower components and is stronger for pi/2 than for the P3/2 orbital. This is in addition 
to the dynamical effect on the upper component due to the Darwin term which is present in 
all the results given here, and that mostly affects the determination of spectroscopic factors 
1^ . It is encouraging that the preliminary data ||22| agree so well with the predictions of 



the fully relativistic calculations and one anticipates being able to make even more stringent 
tests when a finer grid of high-precision data involving other nuclei become available in the 
range 200 < p„ < 400 MeV/c. 
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FIG. 1. Cross section in nb/(MeV sr) (left), Rtl in (middle) and Atl (right) for proton 
knockout from ^^O for the lpi/2 (upper panel) and 1^3/2 (lower panel) orbits, versus missing 
momentum pm in MeV/c. Results shown correspond to a fully relativistic calculation using the 
Coulomb gauge and the current operators: ccl (solid line) and cc2 (dotted line). Also shown are 
the results after projecting the bound and scattered proton wave functions over positive-energy 
states {Jproy, dashed line) and using the asymptotic momenta {Jas, long-dashed line). For the P3/2 
shell a small contribution was taken into account from the nearby 5/2+ and 1/2-1- states known 



from the low-pm data ||26 |. 
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FIG. 2. Rtl (right panels) and Atl asymmetries (left panels) for proton knockout from ^^O 

for the lpi/2 (top panels) and 1^3/2 (bottom panels) orbits. Results shown correspond to a fully 
relativistic calculation using the Coulomb gauge and the current operator cc2 (solid line), a calcu- 
lation performed by projecting the bound and scattered proton wave functions over positive-energy 
states (dashed line) and two nonrelativistic calculations with (long-dashed) and without (dotted) 
the spin-orbit correction term in the charge density operator (see text for details). 
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